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ABSTRACT 

A parallax method to determine transverse velocity in a gravitationally lensed system 
is described. Using the annual motion of the Earth around the Sun allows us to probe 
the local structure of the magnification map that, under certain assumptions, can be 
used to infer the effective transverse velocity. The method is applied to OGLE data 
for QSO2237+0305 and the velocity value is estimated to be about 15 ± lOkms"^ if 
attributed to the lensing galaxy or about 420 ± 300 kms""'^ if attributed to the quasar. 
We find this estimate unreasonably small and conclude that we have not measured a 
parallax effect. We give a short list of properties that a system should possess to allow 
a successful implementation of this method. 

Key words: gravitational lensing — large-scale structure of Universe — galaxies: 
kinematics and dynamics — quasars: individual (QSO2237-I-0305) 



1 INTRODUCTION 

Among the various parameters needed to specify a mi- 
crolensing model of a lensed system one of the most im- 
portant is the overall transverse velocity. This is required in 
order to determine physical parameters from the observed 
temporal properties. Besides gravitational lensing, measur- 
ing the transverse velocities of galaxies is of considerable 
interest for studies of large-sc ale structure and kinematics 
jDekeJIlOgi lBernardeaull2003t) . 

There have been a number of successful attempts to 
use the annual motion of the Earth around the Sun for 
measuring the transverse velocities of microlenses in our 
own Galaxy. The microlensing optical depth remains very 
low within the Local Group and, wherever double lenses 
are not invo lved, can be w ell described by a Schwarzschild 
lens model. iGouldl il992fl developed a convenient formal- 
ism to describe the 'annual parallax' effect in this case, 
that allows one to extract information on physical param- 
eters of the lensing configuration, breaking some of the de- 

J eneracies inher ent to the classical Paczyhski light curve 
PaczvriskilQSd) . The first microlensing event which showed 
a strong annual parallax signal was detected by the MA- 
CHO collaboration in t he first-year data of the group's 
Galactic bulge program jAlcock et alJll995^ and six more 
were found among the longest events over the next seven 
years; two of these are currently among the best stellar 
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mass black hole candidates jBennett et alJl2002^ . Events 
with an annual parallax signal we re also found b y other 
collaborations - E ROS /PLAN ET jAn et alJ l2002h . MOA 
jBond et alJl200lD and OGLE llMadll999l: ISoszvnski et alJ 
l200lHSmith. Mao fc Woznial# 2002T Perhaps the most spec- 
tacular results related to this effect are the detection of 
the first m ulti-peaked microl ensing event by OGLE col- 
laboration jSmith et alJl20o3) . and direct observations of 
a lens based on p arameters determined via this method 
jAlcock et alJl200lD . 

Much effort has also been put into determination of the 
degeneracies present in the annual parallax effect description 
itself. It has been found that constant acceleration of the 
lens or the source can sometime s mimic the parallax signal 
iSmith. Mao fc Paczvriskill2003l) , and there also exists a dis- 
crete degeneracy between jerk and parallax parameters that 
becomes a continuous one when the accu racy of observations 
decreases (lGouldl2004lPark et al .1200411 . These degeneracies 
arise, basically, from the symmetry of the Schwarzschild lens 
that mixes different kinematic effects into a single geometric 
parameter - the distance between the source and the lens 
expressed in the units of the Ein stein-Chwolson radius (e.g., 
'Schneider. Ehl ers fc Falcoll99^ . Some of t hese degeneracie s 
arc broken when the lens is a double object (lAn et al.l2002ll . 
and as lens becomes more complex the above mentioned ra- 
dius loses its unique status in the description of the lensing 
event. The single Schwarzschild lens approximation is, how- 
ever, applicable to most of the microlensing events in the 
Galaxy due to the low optical depth to microlensing in our 
neighbourhood. 

In contrast, microlensing in the images of strongly 
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lensed quasars necessarily takes place in regions of high op- 
tical depth where the structure of the magnification map 
is complex and unknown, effectively rendering it impossible 
to model the light curve in detail. In addition, the intrinsic 
variability that quasars are expected to possess may be an- 
other contaminant of any microlensing signal present in the 
images' light curves. 

For the pa rticular multipl y imaged quasar 
QSO2237-I-0305 (iHuchra et alj [Tosi) in which the mi- 
crolensing phenom enon was observed for the first time 
l|lrwin et al.lll989fl . a number of attempts have been made 
to determine the transverse velocity of the system. The 
small time delay differences between the four images - 
less than a day ijSchneider et aljll98jj) - means that the 
observed fiuctuations, uncorrelated between the images, 
are dominated by microlensing, with negligible intrinsic 
vari ations. 

IWvithe et al.l (Il999ll introduced a method to determine 
the transverse velocity based on the statistics of time deriva- 
tives of the microlensing-induced fiux variations. Although 
their method seems to be the only one that makes it feasi- 
ble to take into account proper motion of the microlenses, a 
number of parameters must be specified in order to make the 
velocity estimate, including the microlensing mass spectrum, 
which, in effect, translates temporal quantities into spatial 
ones. Under the assumptions made in their paper they find 
the transverse velocity to be less than 500kms~^ (at 95% 
confidence level) favouring, depending on the model, values 
in the range 60kms~^ to 400kms~^ - a tighter constraint 
compared to prev iously assumed value of ~ 600kms~^ 
l|Witt fc Madll994l V 

Another approach to estimating the effective velocity is 
to compare the spatial extent of 'quiescent' regions of the 
magnification map models with the temporal extent of the 
periods of steady rise and fall in the actual light curves of the 
quasar images - mostly, image D llGil-Merino et al.ll20oi) . 
However, all these methods are seriously dependent on the 
assumed microlensing parameters and deal with statistical 
properties of the light curves, which are difficult to establish 
with the currently available data. 

In this paper we implement another, rather simple ap- 
proach, underlain by a few natural assumptions. The essence 
of the method is the following: the light curve of an image 
of the lensed quasar is composed of the values of magnifica- 
tion on the observer's path through the magnification map 
(observer's plane). When the region of the map considered 
is small enough and is far from caustic curves, it is natural 
to expand the magnification as a function of the observer's 
position into a Taylor series and restrict ourselves to its lin- 
ear terms. Where this approximation holds, the well-known 
motion of the Earth around the Sun can be used to obtain 
the local values of magnitude and direction of the gradient 
on this magnification map. Combined with the measured 
time derivative, they can be used to estimate the velocity of 
the Sun with respect to magnification map. We apply this 
method to the data available for QSO2237-I-0305. 

In this analysis it is hardly possible to incorporate 
proper motions of the lenses and we will also neglect the 
intrinsic variability of the quasar; therefore it is unreason- 
able to expect a full accounting of all the observed features. 
However, as a (nearly) model-independent estimate of the 
velocity, this approach seems to be an interesting applica- 



tion of the microlensing phenomenon which could be ex- 
tended to different microlensed systems if their properties 
are favourable. 

In the next section we describe the method we use 
for probing the structure of the magnification map and 
obtaining the velocity estimate in greater detail, in sec- 
tion 3 the method is applied to the observational data 
for QSO2237-f0305 obtained by the OGLE collaboration 
( Wozniak et al. ,2000.) and the results are discussed; the 
value for transverse velocity obtained seems to us to be too 
low and we conclude that the method has failed in this case. 
Discussion of desirable properties of lensed systems, which 
may permit successful application of this method, concludes 
this study in section 4. 



2 BASIC APPROACH 

Flux variations in different images of lensed quasars may be 
caused by either intrinsic mechanis ms, such as accretion in- 
stabilities, disc precession etc. (e.g.. |Pringlell99'^ : lKing et alJ 
2004) , or by mechanisms related to the propagation of light 
from the quasar to the observer, which, in the case of lensed 
quasars are dominated by microlensing in the macrolensing 
galaxy. The latter itself could be split into two kinds of vari- 
ations - projected motion of the observer with respect to 
a fixed point in the image plane and 'proper' variations in 
the magnification map itself caused by the proper motions 
of the microlenses. 

In the case of an Earth-bound observer the projected 
motion in turn consists of the constant motion of the So- 
lar system and annual motion of the Earth around the Sun. 
This 'annual parallax' effect in c onnection with microlensing 
of quasars was first considered bv lGrieeer. Kavser fc Refsdall 
(1986) and was shown to be detectable under certain condi- 
tions. However, the authors noted that with the data then 
available it was unlikely that the effect could be detected 
and therefore they did not elaborate on its potential appli- 
cations; for determination of transverse velocities they sug- 
gested usmgsatellite experiments instead (cf. [Refsdal 1963; 
iGoulrf 1995b a). 

The model used in this study is the following: at any 
given time t the observed flux F{t) of a lensed quasar is de- 
termined by the quasar intrinsic flux _fo at the time t — t 
and magnification factor ^ corresponding to the observer 
position (r,t), with gravitational lensing time delay r and 
magnification /i being some functions of the lensing config- 
uration: 

F(t) = Fo{t-r{r,t))x ^i{r,t) (1) 

The observer moves with the Solar system at an unknown 
constant projected velocity v and also rotates around the 
Sun in a well-known annual motion fe(i) as projected onto 
the observer's plane defined as the plane which contains the 
Sun and is orthogonal to the line connecting some fixed 
points in the source and lens planes: 

v{t) =wt + i^{t). (2) 

In the Galactic case, this equation is often written in terms of 
the quantity u(t) - the angular distance between the source 
and the lens divided by the Einstein-Chwolson angular 
radius 6e- 
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(3) 



_ i{t) vt + r,{t) 

U[t) = — — = -:— = TVe z 

oe te 1 au 

The last two equalities define the Einstein radius projec- 
tion onto the image p lane fg and microlensing annual par- 
allax 7r_E = lau/f_B llGouldl f200(i) . The advantage of this 
formulation is that in the case of a Schwarzschild lens the 
magnification depends solely on the magnitude of u, deter- 
mined by the system kinematics, and on Oe, determined by 
the system geometry and lens mass; each of these quanti- 
ties is therefore eissociated with distinct degeneracies in the 
determination of microlensing event parameters. This is not 
the case for the much more complex lenses responsible for 
quasar microlensing, where f_B itself can only be used in cer- 
tain order-of-magnitude estimates. 

If the intrinsic variability can be neglected, time delays 
- which are not independent of the magnification factors - 
drop out of the analysis. This seems to be the case with 
QSO2237-I-0305 as the time delay differen ces between the 
imag es expected in this system are small iSchneider et alj 
Il988l) . while the observed light curves show rich structure 
and do not correlate with each other. Therefore, it is sensible 
to assume that the variability is dominated by microlensing 
and put _fo = const. In this case we can consider a region 
centered at some zero point (ro,to) small enough for the 
Taylor series expansion 

^i{r,t) =^i(fo,^o) + (f-fo)VAi+(^-^o)^+7^ (4) 

to give a negligible remainder, TZ. 

As long as intrinsic variations are not involved it is more 
convenient to deal with astronomical magnitudes. Given the 
times of observations ti and corresponding quasar magni- 
tudes rrii we can attempt to estimate the parameters of the 
representation similar to by fitting it to the data: 

rui - Ami -TZi ^ mo + T{ti — to) + (5) 
X{x,{U) - Xe(io)) + Y{y4U) - y,{to)) 

where Ami is the actual error of the i-th measurement, X 
and Y are the components of the gradient on the observer's 
plane while T is the sum of partial derivative with respect to 
time (caused by microlens proper motions) and 'drag term' 
due to the constant motion of the Solar system: 

T=^+v-Vm. (6) 

All the derivatives correspond to the zero point (ro,to)- 

Since the Earth's motion projected onto the observer's 
plane is an ellipse, the number of parameters fit by © is 
four - two for the straight line corresponding to the aver- 
age motion of the Solar system and two for the amplitude 
and phase of the sinusoid due to the Earth's annual mo- 
tion projected onto the gradient vector of the magnification 
map. This the minimal number of parameters needed within 
this model. We will see in the next section that the quality of 
OGLE data available for QSO2237-I-0305 is such that a four- 
parameter fit fully describes it within observational errors. 
This means that it is not possible to estimate the remainder 
in a proper way and thereby test the validity of the model. 

The major question is whether the partial derivative 
with respect to t in JSJl can be neglected. Although this 
assumption appears to be somewhat customary through- 
out many microlensing studies, random motion of stars in 



the lensing galaxy with respect to each other means that 
it cannot hold pr ecisely and was in fact shown to be an 
oversimplification (iKundic fc Wambsganss"l993^. However, 
the average of this term among different patches of any sin- 
gle light curve should be zero and this can, in accordance 
with common sense, serve as a definition of the bulk velocity 
of microlenses. Therefore, we expect a correlation between 
measured values of T and (X, Y) with coefficients being the 
velocity components: 



T'j ^ ■ Xj +Vy -Yj 



(7) 



Therefore, we need at least three sets of (T, X, Y) to obtain 
a sensible estimate of v. 

The velocity thus obtained represents the sum of the 
transverse velocity in the observer-lens-source system as a 
whole, vo, and the projected local bulk velocity of the mi- 
crolenses at the point where the fc-th macroimage is formed 
with respect to the lens center of mass Ufe transformed by 
the ma crolensing matrix Ak actin g at this point (see (B9, 
BIO) of lKavser. Refsdal fc Staberhl98&) : 



Vfe = Vo -I- AkUk 



(8) 



The Ak matrices are obtained by the macrolensing modeling 
and are usually well-known. 

Macrolensing in the observed systems generally tends 
to magnify images and therefore A matrices effectively con- 
tract the vectors they act upon, while the magnitude of the 
velocities of the microlenses is expected to be smaller than 
the magnitude of the random velocities of the galaxies them- 
selves. Therefore one expects the velocity estimates not 
to vary much from one image image to another and simi- 
larity of these values would provide a strong argument in 
favour of the validity of the model adopted. The arithmetic 
mean can then be used as an estimate for vo which effec- 
tively means requiring ^A^Ufe ~ 0. Except for possible 
differences in weighting, this is equivalent to mixing all the 
images together and fitting one single v in JTJ for all avail- 
able {T,X,Y). 

Whenever Au cannot be considered negligible compared 
to Vo - either due to a chance which exists for the overall ve- 
locity projection vo to be low or because of the high values 
for the individual bulk velocities Uk - additional informa- 
tion on the latter must be invoked, such as the model of the 
expected rotation or other galactic-scale pattern of star mo- 
tion in the lensing galaxy with a subsequent fitting for JTJ 
using Q. 

Another possibility which may be present in this case 
is to utilize the symmetry observed in many gravitationally 
lensed systems which often suggests that the average of bulk 
velocities should be close to zero: ^ Ufe « 0. In this case the 
estimate for the average velocity is 



Vo 



while individual bulk velocities are clearly 
Ufc = Ak^{-Vk ~ Vo) 



(9) 



(10) 



We now apply the above considerations to the gravitation- 
ally lensed quasar QSO2237-)-0305. 
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Hel.JD- 2450000 

Figure 1. Observed light curve of tiie four images of the gravi- 
tationally lensed quasar QSO2237+0305 obtained by OGLE col- 
laboration during OGLE-II phase of the project. Solid lines are 
the fits with parameters given in Table 1 (see text for details). 
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Figure 2. Comparison between calculated T' = Vx ■ Xj -\- Vy -Yj 
and actual Tj values of the coefficient T in joj . The error bars in 
ordinate correspond to uncertainties quoted in Table 1, errors in 
abscissa are the sum of uncertainties in X and Y with coefficients 

^X, Vy. 



3 APPLICATION TO QSO2237+0305 

We use the light curves of four images of QSO2237+0305, 
publi shed by the OGLE collaboration iWozniak et alJ 
l2000l) ^. In the OGLE-II phase of the project over two hun- 
dred measurements were obtained for each of the four im- 
ages of the quasar in more than three years of observations 
(Fig. Ell. The best coverage was obtained over the following 
three periods of time (Hel.JD - 2450000): 650-800, 1250- 
1550 and 1650-1850. In each of these periods the light curve 
of each of the images visually seems to admit representa- 
tion Q, with the exception of A and C images in the sec- 
ond period where distinct peaks, which could be associated 
with a close encounter with a critical line or point, made us 
exclude data with t > 1500 and t < 1380, respectively, from 
the analysis of this period. 

We use the SLALIB package (IWallacd ITooi) to com- 
pute the Earth's position as a function of time and then 
project it on to the plane perpendicular to the line of 
sight toward the quasar, to get and ye- The Ox axis 
is chos en along the COBE dip ole projection onto the same 
plane jLineweaver et al]ll996l) and Oy is such that Ox, Oy 
and the direction to the quasar Oz is a right-hand ba- 
sis (this makes Oxy system apparently left-handed when 
viewed in the image plane). The model is fit to the data 
by the least-squares method - i.e. minimizing the sum 
X^(mo + T{ti - to) + X{xi - xo) + Y(yi - yo) - nnf for 
each of the images in each of the periods independently. 

Results of the fitting are given in table 1. Here we 
present the values of T, X and Y as well as associated er- 
rors. The last two columns show the average fitting quality 
per datum (which square is defined as the sum of residuals 
squared divided by the number of data points) and the aver- 
age observational uncertainty reported by the OGLE group. 
The fits are shown in figure Q as solid lines. 

^ These data are publicly available at 
\http://bulge.princeton.edu/^ogle/ogle2/huchra.html\ 



One can see that the model provides an adequate fit 
for the data given the observational uncertainties. Unfortu- 
nately, this also means that the accuracy of the data does 
not allow for additional tests, and we must admit that a 
simple polynomial in t with the same number of parameters 
can fit the data with almost the same quality. 

We then apply the method presented in the previous 
section to estimate the transverse velocity of the system. 
We first combine twelve sets (T, X, Y) all together regard- 
less of what image they belong to. This is the first option 
mentioned in the previous section valid when transformed 
bulk velocities AhXik are either individually negligible or are 
zero on average. 

To make this procedure reflect the fact that different 
troikas {T, X, Y) have different uncertainties we assigned 
the weight to each of them inversely proportional to the 
square of the average fitting quality from Table 1. In fact 
an unweighted least-squares estimate gives the values of 
Vx , Vy very close (within formal uncertainties) to those of 
the weighted one. 

The best fit obtained is v = (41 ± 7, -1 ± 1.5) kms"^ 
(v = (42 ± 9, -2 ± 2) kms"^ in the unweighted case). Fig.|21 
compares the values of Tj (ordinate) and Tj = Vx ■ Xj + Vy-Yj 
(abscissa). As we have a two-parameter fit here this figure 
should not be viewed as a standard data spread 'cloud'. This 
is rather the best projection of the data points onto one 
single plane determined by v. The fitted values seem to be 
stable to subsampling of the {{T,X,Y)j;j = 1, 12} set. 

The value of the Solar system velocity with respect 
to the background radi ation found by the COBE satellite 
l|Lineweaver et al]ll99(Jl projected onto the observer's plane 
is vcoBE ~ (56.0 ± 0.5,0 ± 0.7)kms~^. Subtracting this 
velocity from the value we obtained we get the transverse 
velocity of the quasar with respect to the lensing galaxy, 
projected onto the observer plane: 'Vq_q = ( — 15 ± 7, —1 ± 
2)kms"^ (vq_g = (-14 ± 9, -2 ± 3)kms"^ in the un- 
weighted case). Taking into account the difference in dis- 



Measuring transverse velocities in gravitationally lensed systems 5 



Table 1. Best fit parameters of the representation |3 for images A - D in each of the three periods (Hel.JD - 2450000): 
1) 650 - 800 (A - D), 2) 1250 - 1500 (A), 1250 - 1550 (B, D), 1380 - 1550 (C) and 3) 1650 - 1850 (A - D). 



T Tn ^ (TO 




± ^mag/yr; 


jC i TY^ cy / CI 1 1 1 


/ wi a <T / 1 1 \ 

1 I lliOj^/ CLU ) 






A 


1 


0.95 ±0.2 


0.057 ±0.04 


0.60 ±0.11 


0.014 


0.014 




2 


-1.41 ±0.04 


-0.098 ±0.01 


0.50 ±0.04 


0.013 


0.010 




3 


-0.37 ±0.07 


-0.073 ±0.01 


0.37 ± 0.09 


0.013 


0.010 


B 


1 


0.91 ±0.4 


0.20 ±0.08 


0.61 ±0.2 


0.027 


0.022 




2 


-0.03 ±0.07 


-0.063 ±0.02 


-0.01 ±0.05 


0.033 


0.032 




3 


0.51 ±0.2 


0.002 ±0.03 


-0.31 ±0.3 


0.036 


0.036 


C 


1 


0.87 ±0.4 


0.22 ±0.08 


0.97 ±0.2 


0.029 


0.035 




2 


3.14 ±0.5 


0.34 ±0.09 


1.66 ±0.3 


0.028 


0.019 




3 


0.43 ±0.2 


-0.026 ±0.02 


-0.33 ±0.2 


0.028 


0.025 


D 


1 


-0.43 ±0.9 


-0.14 ±0.2 


-0.41 ±0.5 


0.068 


0.054 




2 


0.80 ±0.15 


0.035 ±0.04 


-0.29 ±0.1 


0.072 


0.050 




3 


-0.95 ±0.3 


0.01 ±0.03 


0.72 ±0.4 


0.047 


0.047 



tances to the lensinR Ralaxy [zl ~ 0.0394) and quasar [zs ~ 
1.695) iHuchraet alJll985l) . the efTective 'lever arm' ratio is 
about 1:10.5 (in Q, = 0.3, f2 + A = 1 cosmology), and there- 
fore the measured velocity corresponds to either ~ 15 km s~^ 
in the lens plane or ~ 10.5x (l+zg) x 15 kms~^ Ri 420 kms~^ 
in the quasar plane. 

This estimate is curious in both the direction and mag- 
nitude. The direction - which seems to be very close to 
that of the Earth's projected motion with respect to the 
cosmic microwave background - is not, however, a prob- 
lem since the velocity is expected to be aligned with the 
C-D i mage axis based o n an independent quiescence argu- 
ment (I Witt fc Mac fl994'). and the C-D line is almost exactly 
aligned with the COBE dipole projection. Yet the magni- 
tude of the transverse velocity seems to be unreasonably 
small because random individual velocities of galaxies have 
measured dispersion of a few hundreds kilometers pe r sec- 
ond (e.e. iTormen et aljlT99l iBahcall fc Siang"l9 9j) con- 
sistent with theoretical expectations (e.g., Moscardi ni et all 
ll996l:lBahcall. Gramann fc CeJl99'i) thus making a proba- 
bility to obtain a value of a few dozen kilometers per second 
by chance very low (for a Gaussian distribution with one- 
dimensional dispersion of ctc = 300kms~^ the probability 
P of the velocity being less than or equal to v = 15kms~^ 

isP^{;V(2cr|)~io-^). 

This value is also small on the scale of typical ve- 
locities of stars within galaxies which means that the as- 
sumption of Ail to be negligible under which it was ob- 
tained is unlikely to be true. We therefore attempted to 
obtain individual velocities for all four macroimages of the 
quasar. Only three sets {T, X, Y)j - the minimal number 
to make any sensible estimate in this case - were avail- 
able for this task and therefore the uncertainties cannot 
be small. The individual velocities and their formal uncer- 
tainties obtained are va — (61 ± 22,2 ± 4)kms~^, vs = 
(32 ± 20, -4 ± 6) kms-\ vc = (70 ± 110, -7 ± 23) km, s'^ 
and VD = (36 ± 4, -6.9 ± 0.5) kms'^ 

The configuration of the images seen in QSO2237+0305 
shows a high degree of symmetry. We can therefore assume 
that the sum of individual velocities u should be close to 
zero. This suggests to use the second method sketched in 
the previ ous section. To apply it, we use t he macrolensing 
model of lSchmidt. Webster fc Lewis! l)l998h for this lensing 
system to construct microlensing transformation matrices 



Ak for the images. By requiring ^ Uk to be zero, we get the 
value of vo = (48 ± 10, — 2 ± 2) kms~^, while for the random 
velocities we obtain ua ~ (42 ± 20, —4 ± 5) kms~^, ub ~ 
(-49±17, -12±7)kms"\ uc = (20±40, -0.5±20) kms"\ 
u_D = ( — 13 ± 5, 16 ± 2) kms~^; the uncertainties quoted do 
not include the uncertainty of the macrolensing model (i.e., 
those of Ak). 

One can see that both the projected average velocity 
and random bulk velocities at the image locations are very 
low. We are therefore inclined to admit that our method 
described above fails when applied to this system, and the 
parallax signal which must be present in the light curve is 
not detected. Deviations of the light curves from straight 
lines are therefore due to nonlinearity in the magnification 
map and are not caused by the Earth's annual motion — 
that is the TZ term of equation (|1J is not negligible. In the 
next section we describe what sort of systems the parallax 
effect is most likely to be observed in. 



4 WHERE MIGHT THE EFFECT BE FOUND? 

To understand the conditions that would allow one to detect 
the annual parallax effect in multiply imaged quasars, let 
us again rewrite expansion Q about some given point (for 
simplicity, let it denote the coordinate system origin) under 
the assumption that the derivatives with respect to time are 
zero and with the remainder in the Lagrange form: 

m = mo + (vt + fe)Vm+(vt + fe)^I>2 (11) 

Here the first derivative is evaluated at the centre of the 
expansion while D2 represents a sum of second derivatives 
at a point between the centre and vt + f e with certain co- 
efficients of order unity. From this point on the following 
order-of-magnitude reasoning can be adopted. The natural 
scale for the (vt+fe) is the Einstein-Chwolson radius projec- 
tion onto the observer's plane, rs- With lengths expressed in 
these units, the gradient (which equals Vrnffi) is typically 
a quantity of order unity and so is D2f%. 

The strength of the parallax signal is determined by 
TeVm with Ve being an oscillating quantity with the ampli- 
tude of order one astronomical unit and period of one year. 
While fe oscillates, the vt term does not, and therefore the 
ratio of these quantities is optimally measured over a time 
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span of about a year; its maximum value is of order of the 
Earth's orbital velocity over the transverse velocity in the 
system v. The latter is not known for actual gravitation- 
ally lensed systems but it is expected to be about a few 
hundred kilometres per second — that is, much larger than 
the Earth's orbital velocity. Therefore we may now write an 
approximate expression for m: 



m « mo + 
= mo + 



nE[fEVm] + nl[rlD2] 



(12) 



As has already been noted, the quantities in the square 
brackets are supposed to be of order unity, sensible values 
of t are about half a year and v is about a dozen times the 
Earth orbital velocity - (50 — 100) au/yr . The last term in 
this formula represents contamination and the signal is big 
relative to it when the following ratio is large: 



(13) 



(1 - 2)i._iolfE3 



rE3 



That is, the effect is likely to be observed in the systems with 
larger tes = te/ (10^ au) and smaller v^oo = 5/(300 kms^^), 
because the contamination is more rapidly suppresed by in- 
creasing fE (decreasing the microlensing parallax tve) com- 
pared to the useful signal. This is different from microlensing 
by Schwarzschild lenses where no expansion is necessary and 
the second as well as all higher order terms contribute to the 
useful signal. However, te should not be made too large be- 
cause that would mean a very flat magnification map and 
consequently a very weak signal, likely to be lost in the ob- 
servational noise — the only contaminant which is relevant 
to the case of a Schwarzschild lens. Flex points (those where 
the second derivatives vanish) are also highly preferred to 
points of extrema, though this can hardly be controlled. 

The Einstein-Chwolson radius projection f_B = tt^^ au 
onto the observer plane here is t hat of the 'natu ral formalism 
for microlensing' introduced bv iGouldl i2000|) . In a cosmo- 
logical context, it is given by 



4GA/ ^ Dos 



TE = {1 + zl)xI^^Do^ 
which, for a flat Universe amounts to 



(14) 



Te = -1 / (1 + Zl) 



4GM f(zL)f{zs) 
cHo f{zs)-f{zL) 



(15) 



=5.9x10^ {^y\ii+^.)^^mfihy^-- 

(for Ho = 70kms-iMpc-i) with (e.g.. IWeinber3ll972D 



dt 



(16) 



- a monotonic increasing function of z. 

The last ratio in 11511 may be the reason why the effect 
was not detected in QSO2237-I-0305 - close proximity of the 
lens makes the value of f_B relatively low. Placing the lens 
closer to the source, increases the value of fs along with the 
chances to detect the parallax effect but it cannot be placed 



too close as this would reduce the magnification and ability 
to produce multiple images. 

Not much can be said about another player in (I13I I - 
the transverse velocity v. It is generally completely unknown 
and it is for the determination of this quantity that the 
method described here is intended. One of our suggestions 
is to keep away from galaxy clusters, where virialized galax- 
ies are expected to have velocities relatively high compared 
with those of field galaxies. In addition, it should be men- 
tioned that the Solar transverse velocity with respect to the 
CMBR is of the same order as or even greater then the galax- 
ies' velocity dispersion {Bcrnardcau 2003), and therefore the 
expected value of the effective transverse velocity is lower in 
the regions closer to the direction of COBE dipole. 

With regards to the quantity in the square brackets 
in 1131 1 nothing definite is usually known before the observa- 
tions are done. However, there is a way to increase this frac- 
tion by decreasing the value in denominator. Ceteris paribus, 
convolution of the magnification map corresponding to a 
point-like source with an extended profile affects higher or- 
der derivatives of the resulting map in a greater degree than 
lower order ones. Therefore, there is a possibility to suppress 
the contamination from curvature in the intrinsic magnifi- 
cation map relative to the linear signal by observing the 
system in a spectral region where the source looks larger - 
e.g., in broad lines - provided that the smoothing scale is of 
about the same order as the scale of the contaminating sig- 
nal or larger. This effect will mostly be limited by the scale 
at which it will start to suppress the linear term greater than 
allowed by the observational noise. 

Now we can formulate a short list of properties that a 
lensing system should possess to make the detection of the 
parallax effect and determination of the transverse velocity 
possible: 

(i) the lens should be located relatively far away and 
therefore so should be the source 

(ii) the lens and the source should be relatively close to 
each other - though to an extent that still allows multiple 
imaging 

(iii) the lens should not be a virialized member of a mas- 
sive galaxy cluster 

(iv) the system should be located not far from the direc- 
tion of Solar motion relative to CMBR 

( v) either the configuration of the system should be highly 
symmetric or there must be some additional information 
available that could constrain the bulk velocities at image 
locations 

(vi) high degree of symmetry is also an advantage when 
there is a need to constrain the intrinsic variability of the 
source as the time delay differences between images are gen- 
erally low in such systems 

(vii) the lens should be preferably a spiral galaxy with 
a low value of local velocity dispersion at image locations 
because otherwise it would be difficult to link the measured 
velocity of the magnification map with physical velocities of 
microlenses, which, after all, are the actual values of interest 

(viii) where possible, observations in the continuum 
should be supplemented by the observations in broad lines 
or other bands where the source appears larger. 

In addition, standard observational considerations ap- 
ply - the images should be reasonably bright to allow accu- 
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rate photometry while there also should be a way to gather 
some additional information on the lensing galaxy. At the 
moment data of reasonable accura cy and ex tent are only 
available for QSO2237+0305 iWozniak et alJlioOQl) and an 
attempt to determine the transverse velocity in this system 
presented here was not successful. However, as more grav- 
itationally lensed systems are monitored for microlensing, 
some of them may be expected to fulfil the conditions listed 
above, making our method applicable. We suggest that ob- 
servers might usefully focus on such systems. 
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